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Chapter 1 INTRODUCTION 
1.1 General Remarks 


The need to solve increasingly sophisticated soil engineering and soil-structure interaction 
problems has created a demand for better experimental devices and techniques as well as 
better analytical tools, in order to more clear understanding and represent soil behavior. In 
the recent years, considerable experimental and theoretical efforts have been made to 
study and model the soil behaviour under various loading conditions, including cyclic 
loads as caused by earthquakes, traffic and sea wave loads. On the other hand, intensive 
research has also been done on its monotonic behavior, providing abundant valuable 
information on the essence of soil behaviour from very small strains with nearly elastic 
behaviour up to large deformations and flow of the material. 

The rapid development in the recent years of advanced computation procedures for 
nonlinear stress analysis, such as the finite element method, especially emphasized the 
need for appropriate constitutive models, undoubtedly a crucial element in the numerical 
analysis. The development of more advanced constitutive models for soils has also been 
stimulated and thoroughly supported by the increasing supply of high quality 
experimental evidence and a continuous widening of our knowledge of soil behavior. A 
considerable challenge is posed, however, by the fact that soil behaviour is remarkably 
complex, involving various complicated processes and phenomena both under monotonic 
and cyclic loading. 

Constitutive modeling of soils confronts two opposing requirements: on the one hand, the 
model should be able to describe the extreme complexity of soil behavior, and on the 
other hand, it should be comprehensible and simple enough to assure easy 
implementation in numerical analysis and application to practical problems. Whether 
complexity or simplicity will prevail is the first important decision to be made in the 
process of the development of a constitutive model. Yet, this judgment and the choice of 
the relevant aspects of the soil behaviour that should be modeled is probably the most 
critical for the accuracy, effectiveness and applicability of the model. 

It is general trend in elastic-plastic sand models, however, that neither density nor the 
combined influence of density and normal stress is considered as a variable. In fact, a 
sand with different densities is regarded as a different material, and the effects of the 
normal stress are considered to be independent of the density. Roscoe and Poorooshasb 
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(1963) were the first to emphasize the importance of the combined influence of density 
and normal stress on soil behavior. On the basis of both theoretical and experimental 
considerations, Roscoe and Poorooshasb introduced the critical state of the soil as a 
reference for quantification of the effects of the initial state on soil behavior. Here, the 
initial state refers to the void ratio and the mean normal stress of the soil prior to shearing. 
In a similar vein, Been and Jefferies (1985), Bolton (1986), Ishihara (1993) and Verdugo 
(1992) have reported extensive studies on the combined influence of the initial density 
and normal stress on sand behaviour. A distinctive feature of these studies is that 
particular indices and parameters serving as measures for the effects of the initial state on 
sand behaviour have been proposed. The proposed relative Dilatancy index 7 r (Bolton, 
1986), the state parameter (Been and Jefferies, 1985) and the state index I s (Ishihara, 
1993; Verdugo, 1992) integrate the effects of the density and normal stress into a single 
parameter and as such, they are different from the commonly used parameter, the relative 
density D r . Recently, several attempts have been made to use the state parameter and the 
state index I s for sand modelling (Jefferies, 1993; Cubrinivski, 1993; Wood et al.., 1994). 
Jefferies (1993) incorporated the state parameter as a hardening parameter in the Cam 
clay model to describe monotonic behaviour of sand and Wood et al. (1994) used the 
same parameter for modelling drained monotonic stress-strain behaviour associated with 
strain softening. Ishihara and Cubrinovski ( 1998 ) have simulated the undrained 
behaviour of sand using their modified hyperbolic relation along with the stress Dilatancy 
relation based on the energy approach of Roscoe et al 1963, and Gutierrez ( 1993 ) flow 
rule, in plane strain case. 

This study deals with the sandy soils and presents an integrated experimental and 
theoretical investigation on their behaviour under monotonic loading conditions. The 
study contains all the necessary steps in evaluation of a constitutive model for sandy soils, 
including experimental study where the experimental basis for the model is produced, 
then development of a stress-strain relation, formulation of an elastic-plastic constitutive 
model , verification of the model. The present research model follows the framework of 
the state concept, in which a relative initial state with respect to some reference states of 
sand (Critical or Steady state; Quasi-steady state) is used for characterization of sand 
behavior. 
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1.2 On Constitutive Modelling 

Without aiming to review present constitutive models for sandy soils as this point, a brief 
discussion on the general trends in the modeling and basic requirements that a 
constitutive model has to fulfill will be given. To start with, it is important to recognize 
the existence and the need of various types of models since there are considerable 
difference as regards the purpose of the model, the problem that it considers and the 
required level of accuracy and simplicity among others. The development of the advanced 
constitutive models has already been under way for several decades, and these models 
have passed well beyond the speculative or qualitative description of soil behavior. 
During this period, a great number of constitutive models have been proposed involving 
various theories and concepts. Those concepts are either of a micromechanical or 
phenomenological nature. Even tough the principal feature of soils is their particulate 
structure and thus the ultimate adequate representation of the true nature of soils lies in 
the understanding of the microscopic behavior, the complexity of these processes and the 
diversity of the soil field conditions will definitely restrict the quantitative predictions of 
micromechanical models and especially their applicability in the foreseeable future. 
Therefore, it is quite understandable that the majority of the present advanced constitutive 
models rely on phenomenological or macroscopic observations of soil behavior. Only, 
these will be of interest in the present work. 

Recent development in the constitutive modeling of the sandy soils and characterized by 
usage of a variety of approaches and theories including elasticity, viscosity and plasticity, 
as well as combinations of these theories: visco-elasticity, visco-plasticity and elasto- 
plasticity. This trend in the modeling is simply due to the fact that soils exhibit behaviour 
related to these theories. However, as pointed out earlier, it is crucial to select and model 
only the important aspects of the sandy soil behaviour that a constitutive model should 
cover, since that is the only approach which will contribute to the simplicity, 
comprehensibility, and applicability. That this has not been always done in frequent 
constitutive models is clearly indicated by the comment of Scott (1988) on the 
constitutive models presented at the workshop on Constitutive Equations for Granular 
Non- Cohesive Soils (Saada and Bianchini ,1988). 

“Clearly, we cannot go on like this, with the ultimate goal of achieving a thermo¬ 
viscoelastic hypoplastic multiple critical bounding capped yielding endochronic state 
surface models. At some stage, physical intuition has to re-enter the problem....” 
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It has to be recognized , on the other hand, that plasticity models are already proved to be 
capable of capturing the important features of soil behaviour such as Dilatancy, nonlinear 
and irrecoverable deformations, dependency on stress path and inherent and induced 
anisotropy among others. In fact, plasticity theory has been the cornerstone of the 
constitutive modeling of soils for more than two decades. The balance between the 
simplicity and accuracy of a constitutive model seems to be more or less achievable 
within the framework of the incremental theory of plasticity. 

The requirements that a particular model has to fulfill are strongly dependent on the 
purpose of the model and in particular on the problem that it aims to solve. However, in 
general, the model has to satisfy theoretical, experimental and computational 
requirements. In other words, besides the requirements for theoretical rigors, it is essential 
to model the soil behaviour on a physical relevant basis in order to attain a realistic 
representation of the soil behaviour and versatile model. No less important is the fact that 
the number of the model parameters should be reasonable and that their experimental 
evaluation can be easily assessed in conventional tests. Computational efficiency and 
memory requirements are certainly features that can seriously influence the applicability 
of a constitutive model. 


1.3 On Sandy Soil Behavior 

The behaviour of the granular materials can normally be accurately characterized for the 
practical application from the results of the drained tests, since in many problem, the 
material can usually be considered fully drained. However, if the rate of loading is very 
high, such as those resulting from explosions or earthquakes, or the permeability of the 
cohesionless soil is relatively low, as in fine sands or silts, undrained conditions can 
exists. If the material possesses low relative density, significant pore pressure increases, 
may result from these undrained conditions. This in turn may lead to instability, 
liquefaction, and failure earth dam and submarine. Slope failure has been attributed to the 
development of undrained condition in loose, fine sands and silts. These catastrophic 
failures illustrate the importance of understanding the undrained behaviour of sand. 

Undrained monotonic loading of saturated cohesionless materials has been primarily used 
to investigate the mechanics of the soil instability and liquefaction. The requirement for 
soil instability and liquefaction are that a saturated, cohesionless soil exhibits significant 
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contractive volume change tendencies during the undrained shearing with increase in the 
positive excess pore pressure such that continuously decreasing effective stress results as 
shearing progress. This can be achieved by shearing very loose soil at low as well as high 
initial confining pressure, but it can also be realized by shearing the dense specimen at 
high initial confining pressure. When sand is subjected to a monotonic or cyclic shear 
load it exhibits very complex behaviour which is governed by the coupling between the 
shear and the volumetric strains. The link between the distortional and volumetric strains 
is certainly the most distinctive feature of the deformation characteristics of granular 
materials. In accordance with this feature the behaviour of sand is profoundly affected by 
the initial density and fabric of packing, as well as by the initial normal and shear stresses. 
In other words, given sand behaves either in contractive or dilative manner and shows 
different cyclic strength or resistance to liquefaction depending on the initial packing and 
stress conditions. 


1.4 On Triaxial Shearing 

Without any doubt, the triaxial test on cylindrical specimen is presently the most widely 
used procedure to evaluate several important properties of the soil behavior such as: 
strength, dilatancy, stress-strain relationship and anisotropy among others. The triaxial 
test has been used for research as well as for design projects and a very significant part of 
the present knowledge about soil behavior comes from triaxial test results. 

In the conventional triaxial test, the sample is cylindrical in shape, laterally surrounded by 
a rubber membrane and vertically enclosed by rigid end platens through which the axial 
loading is applied. This configuration immediately creates specific boundary conditions 
that, as in any other laboratory test device, are different from the field conditions. This 
shortcoming has been accepted due to the fact that the field conditions can not be exactly 
known nor can they be completed simulated by any laboratory test. In this sense, 
laboratory measurements have always to be seen an approximation of what may actually 
occur in the field. 

Most probably, the world-wide use of the triaxial test in geotechnical engineering is 
because of some of the following advantages of this laboratory test (Baldi et al. 1986): 
relative simplicity of drainage control and measurement of pore pressure, ability to apply 
principal stresses in known directions, ease of measuring axial and volumetric strains, use 
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of solid cylindrical specimens which can be conveniently obtained from standard tube 
samples or easily trimmed from block samples, and versatility of the equipment which 
may be used for a variety of determinations besides triaxial strength and stiffness such as 
for instance, consolidation and permeability parameters, wave velocity, dynamic 
parameters and so forth. On the other hand, the conventional triaxial test has some 
limitations. Perhaps, the two most important are: the principal stresses can not be rotated 
continuously, but only by a sudden rotation of 90°, and it is not possible to have a truly 
generalized stress state, but at least, two principal stresses must to be identical. 
Consequently, the state of stresses is always axially symmetric. 


1.5 Objective and Scope of Present Study 

The main objective of the present study to develop a computational procedure based on 
the state concept to simulate triaxial compressional behaviour of the Bagmati river sand 
subjected to undrained monotonic loading. The intended goal will be obtained through 
these following sub-objectives: 

■ To modify Cubrinovski and Ishihara model based on the state concept for 
torsional simple shear shearing mode to the triaxial sharing mode. 

■ To verify the validity of computational procedure with existing triaxial stress 
strain data. 

■ To calibrate the model for the undrained behaviour of the Bagmati Sand. 

Recent research in soil mechanics is directed towards better understanding of soil 
behaviour and development of advanced constitutive models capable of fulfilling the 
considerable demands related to geotechnical problems. 

In great many circumstances, geotechnical analysis need to accurate modeling of the 
stress strain relation for geologic material. The behavioral modeling of sand requires in 
the bearing capacity and settlement analysis of footing on sand, in liquefaction analysis of 
sand deposits, in sand-pile interaction problems among others. Therefore, a better 
understanding and a more suitable characterization of the undrained response of saturated 
sand are needed. 

Accordingly, the present study is mainly concerned with the simulation of undrained 
response of sand under monotonic loading condition based on state concept. The 
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simulation will be performed for the triaxial compression shearing mode, over a range of 
densities and confining pressure to encompass the fully contractive behavior, partially 
contractive partially dilative behavior, and fully dilative behavior. 

In context of Nepal, there are insufficient soil laboratory facilities and capabilities (for 
limited confining pressure application, and no cyclic loading apparatus). So, there is need 
of simulation researches. 


1.5 Organization of this Thesis 

The thesis is composed of five chapters, the contents of which are briefly summarized in 
the following. 

The introductory part of the study is presented in chapters 1 and 2. Here, chapter 1 
presents background information on constitutive modelling and on the basic concept 
applied in the current constitutive models for sandy soils. In order to better understand, 
the rest of the study, a comprehensive review of the relevant previous studies on sandy 
soils is given in chapter 2. The review includes experimental and theoretical studies on 
sandy soils and places particular emphasis on the essential features of sandy soil 
behaviour. The fundamentals of plasticity theory are briefly summarized as well. 

The theoretical part of the study is presented in chapter 3. Modelling concept and basic 
relations adopted for the present model are elaborated in the chapter. Based on the 
experimental evidence obtained from drained triaxial compression tests, a plastic stress- 
strain relation with density. Stress dependent parameter is developed. The cornerstone of 
the modelling concept is the state concept in the frame work of steady state of 
deformation. The elastic-plastic constitutive model for sandy soils is formulated as well. 
The model is developed in the framework of incremental theory of plasticity and accounts 
for various aspects of sandy soil behaviour for triaxial monotonic loading condition. The 
model is characterized by a single set of parameters for the entire range of densities for a 
given material. 

Chapter 4 covers the verification and validation of the proposed model through 
comprehensive examination in element test simulations. And chapter 5 covers the 
conclusions and recommendations of the present study. 
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Chapter 2 LITERATURE REVIEW 
2.1 Characteristic States of Sand 

There are different characteristic states of sand that have been employed to define an 
index parameter that can incorporate the effect of the initial state on the stress strain 
behaviour of sand. These states include upper reference state and either steady state of 
quasi steady state, and are generally represented by respective lines in e-p’ plane. 

■ Isotropic Consolidation Lines 

For sandy and silty soils, it is conceptually possible to define an upper boundary using 
isotropic consolidation (IC) curve for the loosest state that a given cohesionless soil may 
attain. This curve is equivalent to the normal consolidation curve for clays in the sense 
that defines a limit over which a combination of void ratio and effective mean stress can 
not exist. 

It is important, however, to keep in mind that this curve may be different for different 
sample preparation procedures as shown in fig 2.1.1, Ishihara,1993 ). It is recommended 
that the consolidation curve for the loosest state be evaluated by means of a sample 
preparation by the moist placement (Wet Tempering) method of preparation and then 
saturated. In some soils, it even provides a large void ratio than the maximum void ratio 
evaluated using either the USA or Japanese Codes (Verdugo,1992). 

Similarly to the upper boundary defined by the consolidation line of the loosest state, it is 
possible to define a lower boundary. For the common range of pressure anticipated in 
engineering practice, the isotropic consolidation curve for the densest packing provides 
the lower limit of states in e-p'- plane. For Toyoura sand, both the IC-curves for the 
loosest and densest state are shown in fig 2.1.1 (Ishihara,1993 ). It is readily apparent that 
for the case of the densest state, the volumetric strains caused by isotropic loading are 
very small even though the confining pressure is significantly high. For the practical 
purpose, therefore this curve can be approximated by a horizontal straight line passing 
through the minimum void ratio. This simplification is obviously limited by the pressure 
from which the particle crushing began to be important. According to the test results 
shown in fig 2.1.1, in the case of Toyoura sand, this pressure would be greater than 4 
MPa which includes most of the cases encountered in practice. There are other soils such 
as carbonate containing sand, and decomposed granite sandy soils, however, that are 
much more crushable. For these soil types, the simplification is still valid, but for a small 
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region. Nevertheless, if the particular soil under consolidation seems to be extremely 
crushable, the simplification can be avoided and corresponding isotropic consolidation 
test for densest state performed. 



Fig. 2.1.1 Isotropic Consolidation Curves with different method of sample 
preparation 

The isotropic consolidation curves for the loosest state and the densest states provide a set 
of boundaries for any possible soil state. If attention is drawn to the position of the SSL 
with respect to these boundaries, it is possible to realize that the SSL is indicated very 
close to upper boundary defined by the IC-curve for the loosest state, fig 2.1.2, the area, 
A c , that is associated with those states where soil mass can be potentially unstable small. 
Conversely, fig 2.1.3 shows a material with a SSL closer to the IC-curve of the densest 
state, hence with a dilative area, A j, rather small. 





Fig. 2.1.2 Representation of a soil with small vulnerability against liquefaction 
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Fig.2.1.3 Representation of a soil with large vulnerability against liquefaction 

■ Steady State of Sand 

During the undrained shearing of sand, when it deforms greatly with constant effective 
mean stress, constant shear stress, constant velocity and constant volume, it is noticed that 
the sand attains its steady state for the particular initial state. ( Poulos , 1981 ). During the 
continuous deformation, the original structure of the sand specimen is completely 
destroyed and a new structure "flow structure," which was originally proposed by 
Casagrande ; i.e., due to the process of deformation, the particles become oriented such 
that the shear stress needed to continue deformation eventually reaches a constant value. 
The necessary condition of constant volume and velocity can be achieved in undrained 
shearing performed in strain controlled manner. However, it also occurs in drained tests 
on sand at large strains as well. 

Castro et al.(1976 ); Poulos (1981 ); Been et al. ( 1985 ); Ishihara (1989 ); Sladen and 
Oswell (1989 ); Cubrinovski and Ishihara (2000), have indicated that the steady state 
condition of a particular soil is reasonably unique function of void ratio, fig. 2.1.4. They 
have shown that the steady state is independent of the initial effective confining pressure, 
initial fabric, way of loading, rate of deformation and stress history. However, other 
researchers have pointed out that the undrained steady state strength is not only dependent 
on the void ratio of the soil mass, but is also affected by the effective consolidation stress. 
Konrad ( 1990 )’s interpretation is that there exist upper and lower limits of the steady 
state line. 
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Mean normal stress, p' (kPa) Mean normal stress, bgp' (kPa) 


Presentu(ion of steadj stale tine ia t -p' diagram: (a) e-p’ plot, (b) r-logp' plot 

Fig. 2.1.4 (a,b) Presentation of SSL in e-p’ diagram, Cubrinovski and Ishihara 

( 2000 ) 

When it comes to comparing different soils, however, it has been recognized that both the 
position and slope of the steady state line vary widely in the e-p 1 plot (Castro et al., 1985; 
Poulos et al., 1985; Fear and Robertson, 1995). These studies have shown that the 
variation in the steady state properties is to a great extent due to differences in the grain 
characteristic and grain size composition of soils. The vertical position of the steady state 
line is chiefly influenced by the grain size distribution, whereas the slope of the line is 
affected by the grain shape (Poulos et al., 1985). It has also been recognized that fines 
have some influence on both position and slope of the steady state line (Been and 
Jefferies, 1985; Sladen et al., 1985; Pitman et al., 1994). 

The investigated steady state lines have been obtained from triaxial compression tests on 
reconstituted samples of clean sands and sands with a non-plastic fines content of less 
than 20%. 

The assumption that the steady state strength is a sole function of void ratio leads to the 
conclusion that such line should be unique for given sand irrespective of the method of 
testing and sample preparation method. But several experimental results show that it is 
not always true and pose serious question about its uniqueness. Several researchers in 
detail have studied different factors that can affect the apparent position of the steady 
state line in case of clean sand. Uniqueness of steady state line seems to be studied mainly 
with respect to the factors such as 
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• Mode of deformations 


• Method of sample preparation 

• Drainage conditions 

• Forms of consolidation 

• Bedding orientation etc., 

In literature contradictory conclusions regarding uniqueness of the steady state line with 
respect to mode of deformation can be found. Results obtained by Ishihara et al (1998) on 
sand fraction of Masado soil, Cubrinovski et al (1999) on Kasumigaura sand reveal that 
mode of deformation has indeed no effect on the position of steady state line. 




Fig. 2.1.5 Effect of mode of deformation on the position of SS line on sample 
prepared by a) Wet Tamping b) Water Sedimentation, Ishihara et al. (1998) 

Fig. 2.1.5a, 2.1.5b show steady state line obtained in triaxial compression and extension 
from the sample of Masado soil prepared by the method of water sedimentation and wet 
tamping or moist placement. 

These test results strongly support the fact that steady state line is unique with respect to 
mode of deformation. On the other hand results obtained by Ishihara and Yoshimine 
(1998) on Toyoura sand, Reimer and Seed (1997) on Monterey sand show that mode of 
deformation has distinct effect on its position. Steady state lines obtained in triaxial 
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compression, extension and simple shear in case of Monterey sand are displayed in figure 
2 . 1 . 6 . 



Fig 2.1.6 Non-uniqueness of SS line with respect to mode of deformation in case 
of Monterey Sand 


From this figure it is clear that undrained strengths measured in triaxial compression in 
considerably higher than that obtained in triaxial extension. Simple shear seems to show 
even lower value than either triaxial test. In this study, only moist placement method of 
sample preparation was employed. On the contrary simple shear showed intermediate 
value higher than triaxial extension and lower than triaxial compression in case of 
Toyoura sand prepared by dry deposition method, Yoshimina and Ishihara, 1998. 
Negussy and Islam, 1994 also reported that the steady state strength of tailing sand 
measured in triaxial compression mode are different and higher form extension mode. 
Thus in the case where different steady state line is observed higher steady state strength 
always seems to be associated with compression mode and the difference in the value 
obtained in these modes increases as the relative density of the sand decreases. It is, 
therefore apparent that application of steady state line obtained form triaxial compression 
test would lead to an inappropriate assessment of liquefaction potential (figure 2.1.7 ). 
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Fig. 2.1.7 Undrained behaviourof Toyoura sand in triaxial compression, extension 
and simple shear tests ( Yoshimine et al., 1999b) 

The second factor considered here is the effect of sample preparation method or mode of 
deposition. Steady state line obtained from the reconstituted sample of Masado soil 
prepared by water sedimentation and wet temping is shown in figure 2.1.8a and 2.1.8b 
respectively. It is evident from the figure that similar strength can be obtained regardless 
of the method of sample preparation employed in triaxial compression and extension 
modes. 
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Previously this fact was revealed by Verdugo (1992) where three methods of sample 
preparation was applied and unique steady state line of Toyoura sand was obtained in 
triaxial compression. However, there exist some experimental results that indicate effect 
of fabric on its position, DeGregorio (1990). 

Steady state line appears to be not affected by drainage condition applied during shearing. 
This conclusion is based on the experimental results obtained on Toyoura (Verdugo, 
1991) and Monterey (Reimer and Seed, 1997) sand prepared by moist placement in 
triaxial compression. It is not known whether similar conclusion could be drawn in case 
of other modes of deformation and deposition. However experimental observation in case 
of sand containing silt revealed that the steady state line obtained from drained and 
undrained test coincides high at higher pressure but diverges at lower pressure, 
Yamamuru and Lade, 1998. It is to be noted that the strength of sand at steady state 
decreases with increasing fines content. Similarly form of consolidation; over 
consolidation ratio and bedding orientation have practically no effect on the apparent 
position of the steady state line, Kato et al (2001) and Verdugo (1992). From the above 
discussions, it can be concluded that the issue of uniqueness of steady state line is 
controversial especially with respect to modes of deformation employed. So far, the 
existing literature fails to explain under which condition steady state line is unique for 
given sand and hence merits further investigation. Although above discussion put doubt 
on the uniqueness of the steady state line, yet it is assumed that factors considered above 
have no influence on its position in order to check its appropriateness in defining flow 
and non-flow conditions of sand and in modeling the behaviour of sand. 

■ Quasi steady state 

When a loose sample of sand is subjected to relatively high confining stress and sheared 
in monotonic undrained manner it exhibits the contractive behaviour at the beginning but 
tends to exhibit dilative behaviour at the later stage of shearing. The point at which the 
behaviour changes from contractive to dilative has been defined as the phase 
transformation (Ishihara et al., 1975). During the shearing action, a temporary drop in 
shear stress may or may not occur. If the temporary drop in shear stress takes place at the 
point of phase transform over a limited range of shear strain then the sample is said to be 
in quasi steady state (Alarcon-Guzman,1988 )(Been et al.,1991 ). Hence, the quasi steady 
state can be regard as a particular case of phase transformation where the transition from 
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contractive to dilative behaviour takes place with a drop in shear stress. The occurrence or 
non- occurrence of the quasi steady state primarily depends on the initial confining stress 
applied prior to shearing. The strength of the sand at this state is considerably lower then 
at the steady state, fig. 2.1.13. However, the friction angle mobilized at SS and QSS has 
been considered to be a unique material property and does not depend on the initial state, 
direction of principal stress and mode of loading; Vaid et al., 1990, Ishihara, 1993, 
Uthaykumar and Vaid, 1998. 



Fig 2.1.10 Three characteristic lines in the e-p’ plane 

Unlike the SSL, the position of the QSSL is less arguable. It indeed depends on the 
method of sample preparation, Verdugo, 1993, Zlatovic and Ishihara, 1997,.This state of 
sand is inherent in the fabric of sand as different methods of sample preparation accounts 
for different fabric. In the references cited, it can be seen that QSS has been studied 
mainly in triaxial compression mode. Like SSL, this is also unaffected by the factors such 
as form of consolidation, over consolidation ratio etc. Verdugo,1992 and Kato et al.2001 
respectively. 
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■ Upper Reference State 

Experimentally, it has been observed that given sand at same particular value of void ratio 
starts to exhibits zero residual strength upon undrained shearing. The void ratio at which 
zero residual strength upon undrained shearing first occurs is called limiting void ratio or 
threshold void ratio. Any other values of void ratio higher that threshold void ratio is 
always associated with zero residual strength when sheared in undrained condition. 
Theoretically, this value of void ratio is determined by extending either SS or QSS-line 
up to the void ratio axis in the e-p’ plane. Intersection point at the void ratio axis gives the 
value of threshold void ratio, Ishihara,1993 and Verdugo,1992. It can be seen, in general, 
that SS and QSS line are linear in the e-p’ plane, when the effective stress is greater than 
10 Kpa. Below this value, these lines are characterized by steep slopes. Hence, the 
extension of the SS or QSS line up to the void ratio axis requires adequate attention. A 
line passing through the threshold void ratio, e Q and parallel to the effective stress axis 
has been termed as upper reference line ( UR-Line ). When this line intersects the IC-line 
for the loosest state at some particular value of effective stress then, beyond the 
intersection point, IC-line for the loosest state serves the purpose of the UR-line. 

2.2 Index Parameters 

■ State Parameter 

In particular, confining pressure modifies material behaviour of sands to the point that 
even dense sand, if tested at sufficient high confining pressure, will behave similarly to 
loose sand. Therefore, properties of sand can not be expressed in term of relative density 
alone; a description of stress level must also be included. 

When same sand is tested under very different combination of void ratio and mean 
effective, behaves similarly, if test conditions assure an equal initial proximity to the 
steady state. Been and Jefferies (1985), defined the state parameter, which is defined as 
the void ratio difference between the initial sand state and the steady state condition at the 
same effective stress. 

y/ = e-e Q (2.2.1) 

A single parameter, the state parameter, can be used to describe much of the behaviour of 
the granular material over a wide range of stresses and densities. While this does not 
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eliminate the influence of either density or confining pressure on the behaviour of sand, it 
properly places emphasis on the fact that it is a combination of these parameters that is 
physically relevant to the description of the granular materials. The parameter can be used 
to normalize large strain behaviors where the influence of the initial fabric is small. 
Anisotropic fabric is regarded as a second parameter which, in combination with the state 
parameter, will allow a reasonably complete description of sand behaviour at both small 
and large strains. The physical consideration for an appropriate state parameter for sand 
are therefore what is the appropriate combination of the void ratio and stress and what is 
the appropriate reference condition. The selection of a reference condition on which to 
define state requires that the reference condition should have a unique structure which is 
not influenced by the original test conditions. This can be achieved only at steady state. 

■ State Index 

The State Index, I s , proposed by Ishihara ( 1993 ) and Verdugo ( 1992 ), is used to 

quantify the link between the relative initial state and the normalized stress-strain curve. 
Two characteristic states of sand are employed in the definition of Is, i.e. QSS/SS-line 
and UR-line. The UR-line is associated either with the threshold void ratio, e Q , above 
which the initial states lead to a zero strength, or with the isotropic void ratio e and mean 
normal stress p, the state index is defined as, 



Where, eg is the void ratio of the QSS/SS-line and ey is the void ratio of the UR-line, 

both at the initial mean normal stress p. Thus, the state index is a direct measure for the 
relative initial e-p state with respect to the QSS/SS-line and UR-line. As shown in 
Fig.2.2.1, when the initial mean normal stress is p < p Q ,the upper reference line is a 
horizontal line corresponding to void ratio e 0 . On the other hand, for initial mean normal 
stresses in excess of p 0 (p>p 0 ), the isotropic consolidation line for the loosest state (IC- 
line) is used as the upper reference line. The void ratio e 0 is defined as a threshold value 
at or above which the initial states are associated with a zero steady state strength. In 
other words, upon monotonic undrained shearing of samples with void ratios e > e 0 , both 
shear and normal stresses drop to zero and the soil exhibits ‘static liquefaction’. The void 
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ratio e 0 can be evaluated as the point of intersection between the QSS -line and y-axis. 
Alternatively, the steady state line, since for any initial state with e > e 0 the steady state 
and the quasi steady state are coincident. The mean normal stress p 0 denotes the normal 
stress at which the e 0 -horizontal line intersects the isotropic consolidation line for the 
loosest state. The above expression of I s yields the conditions I s = 0 and I s = 1, for the 
initial states at the UR-line (e = e v ) and at the QSS/SS-line (e = e Q ), respectively. The 
choice of the reference states in the definition of I s is elaborated by Ishihara (1993). It is 

important to note that the above frames of reference are inherent in the fabric of the sand 
formed during its deposition, and therefore, they are fabric-dependent. 



Fig. 2.2.1 Definition of State Index (Ishihara, 1993 ) 

It is seen that initial states that have similar I s -value have similar stress strain curves, as 
the I s decreases, both the initial shear modulus and peak stress ratio decrease as well; the 
degradation of the stress strain curve ceases once I s approaches the UR-lines( I s =0 ). It 
is also seen that a rational relationship between the stress strain parameters and state 
index can be established when physically meaningful are used. The strong correlation 
between the strain stress parameters and I s is enhanced by the fact that employed 
characteristic states are not mere reference lines but rather principal states of sand 
behavior. Employing the state index as a current variables enables us to gradually change 
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the characteristic of the stress strain curve according to the current e-p state and achieve a 
density-dependent influence of the mean normal stress. 


■ State Pressure Index 

Dafalias( 2002 ) has defined the state pressure index, 1 p as the ratio of the current mean 
pressure to the mean pressure at steady state that correspond to the current void ratio. 


/ 


P 


P 

PQ 


(2.2.3 ) 


The concept of state pressure index is introduced as an appropriate state variable for the 
constitutive modeling of sand. The mean pressure at critical state, pg , is determined 

from the current void ratio, e using the critical-state line. For example, if the-critical state 
line is given by 

e Q = r - A • log(pg) (2.2.4) 

then, substituting e for eg , we have 

r—e 

p Q = 10 * ( 2.2.5 ) 

The definition of I p is valid for both drained and drained behavior. For undrained 
conditions, because e is constant, pg is constant. The relationship between different 
state parameters can also be found. In general, if the critical-state line is given as pg = f 
( eg ), then the current pressure p can be related to a critical void ratio eg , and the critical 
pressure p c can be related to the current void ratio e. Thus I p = p / pg = f ( eg) / f (e). 

The state pressure index, I p , can be related to the state index, I s as: 

I P = — = I s ~e ( 2 - 2 . 6 ) 

PQ 
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2.3 Undrained Monotonic Behaviour of Sand 


The ‘Undrained’ soils are two-phase particulate materials in which the voids between the 
particles are full of water. In addition, the permeability of the material may be sufficiently 
low or the loads applied so quickly that pore water pressures that are generated have no 
time to dissipate during the time to dissipate during the time-scale of the analysis. When 
saturated sand with sufficiently lower value of void ratio is sheared in a monotonic 
undrained manner, it exhibits a peak value of shear strength at the beginning. As the 
deformation continues, a decline in the shear strength takes place due to the generation of 
increasing pore water pressure. This phenomenon can be well described with reference to 
Fig.2.3.lb where characteristic stress strain curves are respectively shown. 


(a) 
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(Strain hardening) 
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Fig. 2.3.1 Undrained behaviour of sand (a) Effective stress path, (b) Stress strain 


curve 


The response in such condition can be explained in three different ways. Following peak 
shear strength, very loose sand loses its strength and exhibits strain- softening behavior. 
This strain softening behaviour continues until a constant value of pore water pressure is 
reached. The response associated with this type has been referred to as true liquefaction 
or flow deformation. In second type, sand with some higher value of void ratio shows 
strain softening followed by strain-hardening response. The strain hardening takes place 
as a result of the generation of negative pore water pressure during which sand regains 
some strength and restores stability. Again this hardening phase continues until constant 
value of pore water pressure is reached. One salient feature of this response is that it 
shows two different states viz. quasi steady state and steady state at moderately large 
shear strain and at large shear strain respectively. 

Due to the exhibition of strain softening as well as strain hardening behavior, this type of 
response is generally knows as limited liquefaction or limited flow. Third type is 
associated with medium dense to dense sand in which only strain hardening is generally 
observed and ever increasing shear stress is required to produce shear deformation. On 
further shearing, steady state is reached where no further increment in effective mean 
stress as well as shear stress takes place. 

It is clear that steady state is generally achieved at sufficiently large deformation 
depending upon the relative density of sand. The shear strength mobilized at this state is 
referred to as strength or residual strength of sand. However, in some cases, the strength 
mobilized at quasi steady state is considered as the strength of sand. 

The several phases of the stress strain curve in case of medium loose sand are shown, in 
general, up to the range of 0.1% , stress strain curve is linear and elastic. Peak value of 
shear strength occurs when the deformation is in the range of 0.1-2%. Thereafter strain 
softens and reaches a minimum value. The occurrence of this minimum value is observed 
in the moderately large range of strains (2-15%). Finally strain hardens and reaches the 
steady state when the deformation is in excess of 15% 

Hence complete change in the original structure of the sand specimen is regarded as the 
sufficient criterion for sand to be in steady state deformation the necessary condition of 
constant volume and velocity can be achieved in undrained shearing performed in strain 
controlled manner. However it also occurs in drain test on sand at large stain as well. 
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Under this state, the condition of constant shear stress and constant normal stresses are 
also satisfied to large extent. This ultimate condition corresponding to the actual steady 
state and is always achieved, irrespective of whether the specimen is contractive or 
dilative. 


2.4 Factors Influencing Sandy Soil Behaviour 

The behaviour of sandy soils is affected by a number of factors including density, 
effective stress, shear stress, strain history, initial fabric, and orientation of principal stress 
directions and relative magnitude of intermediate principal stress among others. A great 
many experimental studies have been conducted on sandy soils in order to study the 
effects of these factors and to quantify their influence on soil behaviour. To demonstrate 
the influence of the most important factors and to illustrate the essence of the behaviour 
of sandy soils, some typical data obtained from studies conducted on cohesionless soils 
under drained conditions will be presented. The undrained characteristics of the 
cohesionless soils will be discussed in the subsequent section. 


■ Influence of Initial Void Ratio 

From the early studies in soil mechanics it has been recognized that the initial density or 
void ratio and effective stress have most important influence on the stress-strain an 
volume change characteristics of soils. Thus, traditionally in soil mechanics the general 
patterns of soil behaviour have been related to parameters associated with density of the 
soil, such as relative density and void ratio. Characteristics of behaviour of cohesionless 
soils associated with different initial densities, as those typical for loose sand and dense 
sand, are now well established and supported by abundant conclusive experimental 
evidence. Herein, influence of the initial density on behaviour of sand will be illustrated 
with the data presented by Koerner(1970), obtained from conventional drained triaxial 
tests. 

There is significant influence of the relative density on both stress-strain and volume 
change behaviour. The stress-strain curves are nonlinear, with gradual change of the 
stiffness and peak stress ratio depending on the relative density of the soil. Dense soil 
exhibits a peak in the stress-strain curve and subsequent strain softening post-peak 
behaviour, while the loose soil is characterized by continuous increase in the stress ratio 
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until the ultimate stress state is reached. However, a unique residual stress, independent of 
the relative density, is reached at a certain axial strain. The stiffness of the stress-strain 
relationship is greater (or the slope is steeper) with increasing relative density. 

The volume change behaviour in these tests is characterized with overall dilative 
response. There is gradual increase in the compressibility as the initial relative density 
decreases, with the dense soil being slightly compressive in the initial part of the shearing 
and strongly dilative as the shearing proceeds. The loose sand, on the other hand, exhibits 
strong compression and relatively poor dilation. 

■ Influence of Mean Effective Stress 

Mean effective stress is another major factor with significant influence on behaviour of 
sandy soils. It has been subject of numerous studies on cohesionless soils utilizing drained 
conventional triaxial compression tests (Lee and Seed,1967; Vesic and Clough, 1968; 
Fukushima and Tatsuoka, 1984 and Kolymbas and Wu, 1990 among others), triaxial tests 
on large, squat samples (Hettler and Vardoulakis, 1984), plane strain compression tests 
(Tatsuoka et ah, 1986a) and torsional shear tests on hollow cylindrical samples(Tatsuoka 
et al., 1986b) among others. Except for the torsional shear tests on hollow cylindrical 
samples, other conventional tests mentioned above are performed at constant effective 
confining stress oc, while the mean effective stress is not constant but rather increasing 
with the deviator stress. Triaxial tests under constant mean effective stress, e.g. as those 
reported by Kondner and Zelasko (1963) and Vesic and Clough (1968), have seldom been 
performed. In general, however, there is conclusive evidence on the influence of the mean 
effective stress, although some details remain to be clarified. 

Lee and Seed concluded that an increasing confining stress had three effects on the 
behavior: it reduces the brittleness of the stress - strain curve, it increases the strain to 
failure or to peak stress ratio, and it decreases the tendency to dilate. They further 
concluded that the stress - strain characteristics of dense sand at high confining pressure 
are not unlike those of loose sand at low pressures. 

One of the important conclusions in the study of Lee and Seed is that at high confining 
pressures, the effects of particle crushing on behavior of cohesion less soils is significant, 
and therefore, they must be properly considered. 
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2.5 Elastic Behaviour of Sandy Soil 

Experimental evidence on the behaviour of cohesion less soils has sown that practically 
both elastic and plastic deformations occur from the beginning of loading. A similar 
conclusion has been derived by Mindlin from the theoretical consideration of two elastic 
spheres pressed together, which states that slippage must occur even for the smallest 
tangential force (Hardin, 1978; Tatsuoka and Shibuya, 1991). Thus, elastic and plastic 
deformations are in face always in parallel, and therefore, it is difficult to isolate them 
experimentally. Even though calculation of the elastic strains is related to considerable 
experimental difficulties, and their magnitude is relatively small in the presence of plastic 
strains, their precise evaluation is still extremely important. Basically, there are three 
reasons for such a need: first of all, any accurate calculation of the behaviour of sandy 
soils at very small strains (e.g. less than 10" 5 ) relies on the elastic parameters since at such 
small strains the elastic behaviour dominates; on the other hand, in the formulation of the 
flow theory of plasticity only the plastic strains are involved, and therefore, it is necessary 
to deduct the elastic contributions from the experimentally obtained strains; finally, 
constitutive models must include elastic strain components in order to invert the stress- 
strain relations. Hence, the elastic behaviour has to be appropriately represented in the 
constitutive model in order the assure accurate simulation of the soil behaviour over the 
whole required range of strains, from small strains, up to large strains and failure of soil. 

Despite the difficulties related to the determination of the elastic parameters there are 
many experimental procedures for their evaluation, including laboratory tests, field 
loading tests and field shear wave velocity measurements. As regards the laboratory tests, 
it is considered that at least a good estimate of the elastic moduli can be obtained from the 
following: the initial slope of the virgin stress-strain relation (Tatsuola and Shibuya, 
1991); the initial slope of the unloading curve (Cole,1967; Lade and Nelson,1987); from 
the slope of the unloading-reloading curve (Duncan and Chang,1970); and from 
proportional loading tests (constant stress ratio) with decreasing effective stress and shear 
stress (Rowe, 1971) among others. All these methods for determination of the elastic 
parameters are based on the assumption the soil exhibits purely elastic behaviour or that 
at least this behaviour is dominant under those loading conditions. With respect to the 
above mentioned methods, it appears that the first two can provide much better 
information about the elastic behaviour of sandy soils, provided that both stresses and 
strains are accurately measured at very small strain level, while the two later methods 
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provide relatively rough data for the elastic properties since these loading conditions 
involve plastic strains of comparable magnitude to the elastic ones. 

Based on the experimental results Tatsuala and Shibuya (1991) concluded that the 
behaviour at very small strains is linear and recoverable and they set a ‘limit-strain’ at 10" 
5 below which the behaviour is considered as linear and elastic. Another important 
observation by the same authors is that the elastic region does not expand with increase of 
the applied maximum stress level, and therefore, they suggested using ‘kinematic 
hardening modeling with fixed and rather small size of elastic region’ in order to 
realistically represent the behaviour of sands. The shear modulus at unloading is equal to 
the shear modulus measured at the initial part of the curve only at the very beginning of 
the unloading part, thus demonstrating that the elastic region is very small. 

It is generally accepted that for practical purposes the elastic behaviour of sandy soils can 
be considered as isotropic, and that all anisotropic behaviour should be referenced to the 
plastic part of the response. This assumption simplifies the matters considerably and 
requires only two independent parameters for describing the elastic behaviour, and these 
may be Young’s modulus E and Poisson’s ratio v or the shear modulus G and the bulk 
modulus B. 

Experimental results indicate that elastic properties are dependent on the state of the soil, 
i.e. they may be expressed in terms of void ratio or density and mean effective stress. 
Thus, it is customary to assume that Poisson’s ratio v is constant, whereas the Young’s 
modulus E, the shear modulus G and the bulk modulus B are represented as functions of 
void ratio and density. It should be emphasized that the elastic moduli are related to the 
density of the soil, or to void ratio e, rather than to the relative density D r , as pointed out 
by Iwasaki and Tatsuoka (1977) and Alarcon-Guzman et al.(1989). 

Most of the expressions for elastic shear modulus G presented in the literature have the 
following general form: 

G = AF{e)p n (2.5.1) 

Where A is a constant, F(e) is a function describing the influence of void ratio, p’ is the 
mean effective stress and n is an exponent which is usually taken as 0.5. Following this 
general form, Iwasaki and Tatsuoka (1977) and Iwasaki et al. (1978) proposed a set of 
equations for sandy soils where G is expressed as: 
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G = A(y)B- 


(2.5.2) 


(2.17 e) n (/) 
l + e 

In which A(y) and n(y) vary from 900 to 700, and 0.4 to 0.5, respectively, for the strain 
range between 10 ~ 6 to 10~ 4 . The variation of these parameters with strain indicates that 
even at these small strains, the behaviour was yet not purely elastic. However, for 
practical purposes the values proposed for the smallest strain at 10" 6 can be assumed to 
represent the elastic behaviour. As regards the parameter B, Iwasaki and Tatsuoka found 
that for uniform clean sands, B is equal to unity and independent of the grain size in the 
range between D 50 =0.16-3.2 mm. On the other hand, they concluded that B is strongly 
affected by the non-uniformity of the sand and by the fines content such that B decreases 
as sand is more non-uniform and as the presence of fines increases. 

2.6 Fundamentals of Plasticity 

■ Basic assumption 

The theory of plasticity has already been established as the cornerstone of the constitutive 
modeling of soils, and has been accepted as a frame work within which soil behavior can 
be described providing the balance between the simplicity and accuracy of a constitutive 
model. 

It is assumed that the increment of total strain df ,, can be divided into its elastic d; y e 
and plastic d£ij P components. Mathematically, 

d ?ij = de i/ + de / ( 2.6.1 ) 

Essentially , the plastic strains remain after the removal of stresses while the elastic strain 
are recoverable. The elastic strains and the plastic strains are calculated separately, by an 
elastic stress-strain law and by the theory of plasticity including yield surface, hardening 
rule and flow rule as its constitutive elements. 

■ Yield Criteria 

A yield surface defines the limit of elastic behavior for a material under a combination of 
stresses. Inside this surface, the behavior is assumed to be elastic and therefore , strains 
can be recorded upon the removal of the stresses. The second equally important role of 
the yield surface is to denote the stress space where plastic strains occurs. 
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Based on experimental evidence on soil behavior , various yield surface have been 
proposed , having different size , shape and rules for their evolution. The equ. of yield 
surface, 

/ =/ ( o’y ) ( 2.6.2 ) 

■ Flow Rule 

The relationship between elastic strain increment and stress increment is expressed by 
Hook’s law. For the yielded material satisfying progress consistency condition, a flow rule 
relates the plastic strain increment to the current stress. 

Besides that, it also specified the direction of the plastic strain increment in terms of the 
plastic potential g, 

dp 

de ij e = <L> (2.6.3) 

a 

where , L is a quantity called loading index. It determines the magnitude of plastic strain 
increment and depends on current state of stress and loading history. The bracket 
enclosing L is called Macaulay’s bracket and has the following meaning , 


<L> 

- 0 

for 

L 

< 

0 

<L> 

= L 

for 

L 

> 

0 


Flow is said to be associate if f = g otherwise, it is considered as non-associative. In fact, 
the purpose of plastic potential is to determine the relative magnitude of the plastic strain 
increment. However, it is generally assumed that the plastic potential is a function of 
current stress level alone and is neither affected by the loading history nor by the 
magnitude or direction of the stress increment. 

■ Hardening Rule 

While yield criterion and flow rule specify the elastic limit and direction and magnitude 
of plastic strain increment respectively, a hardening rule is generally introduced to specify 
the size, position and shape of the yield surface. In general, the size position and shape of 
the yield surface are not fixed and may change in the stress space during loading. It is 
necessary that, as stated earlier, the stress point must always lies on the yield surface and 
hence the size position or shape of the yield surface are required to be determined. For 
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this purpose, a hardening rule is introduced. In practice three types of hardening laws 
exist. 

■ Isotropic Rule 

In this the type of hardening law, the yield surface expands uniformly everywhere about 
an axis. Isotropic hardening takes place when the initial yield surface in the stress space 
expands uniformly without changing its shape during the plastic flow and the center of 
the yield surface remains the same during the expansion. In the other word, yield surface 
expand without distortion and translation during the loading process as shown in fig. 
2.6.1. In general, the size of the yield surface can be represented by a single parameter 
such as k, the hardening parameter. Whenever, the stress point touches the yield surface 
and tries to move beyond it, the yield surface expands with the amount of the expression 
given by dk as determined from consistency conditions. The yield surface then grows 
uniformly until it touches the failure surface. The mathematical expression for the yield 
surface can therefore be written in the form 

f(<T ij ,k)=f(<T ij )-k{k) = 0 (2.5.6) 

CTi 



Fig. 2.6.1 Concept of Isotropic Hardening 

Any parameter that is increasing monotonically can be used as isotropic hardening 
parameter. Plastic volumetric strain ( Schofield and Worth, 1968 ) and plastic shear strain 
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( Poorooshasb and Pietruszcak , 1985 ) are the most common in use. In the literature, a 
linear combination of plastic volumetric and shear strain has been used as a hardening 
parameter (Nova and Wood, 1979). Lade 1975 used plastic work as a hardening 
parameter in his work. Isotropic hardening rule is simple but can not fully represent the 
response of the soils under complex loadings. 

■ Kinematic Hardening 

In this type of hardening rule, it is generally assumed that the loading surface or yield 
surface only translates as rigid body in the stress space during plastic deformation. It 
means that the size, shape and orientation of the initial yield surface are maintained , if 
the yield surface are not allowed to rotate . It is the kinematic hardening rule that provides 
the simple means of accounting for the Bauschingers effect as shown in fig 

CTi 
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Fig. 2.6.2 Concept of Kinematic Hardening 

The yield surface for the kinematic hardening material is generally expressed by 
/(try , «ij) = f (CTij - oty ) - AT-0 (2.5.7) 

where k is a constant. a ,, is called a back stress and gives the coordinates of the center of 
the surface. The hardening response is represented only when this back stress varies 
during plastic deformation. Hence, key step in this rule is the determination of back 
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stress. For the sake of convenience, a linear dependence of d(X y on d< ,j P can be assumed. 
This is known a sPragers hardening rule and can be expressed as 

rfoij = ci d£ ij p ( 2.5.8 ) 

where cj is a material constant and may be a function of state variables such as k. 

• Mixed Hardening Rule 

A combined isotropic and kinematic hardening rule can be used in order to model the 
response of materials to cyclic loading where progressive decrease in the magnitude of 
plastic strain occurs during successive loading. This combination results into a more 
general form called mixed hardening rule. Obviously , if this type of hardening rule is 
assumed both expansion and translation of the yield surface takes place. Mathematically , 
such type of hardening characteristics can be expressed by the equation of the form 

/"(cry , fXij, k ) = /(cry , ay ) - K (k) = 0 ( 2.5.9 ) 


&1 



Fig. 2.6.3 Concept of Mixed Hardening 

In this case the expression of the loading surface is measured by K (k) whereas translation 
by ay . But it still retains its original shape as shown in fig. 2.6.3.With the mixed 
hardening rule , different degree of Bauschinger's effect can be simulated by adjusting 
two parameter , K (, k ) and (Ty . 
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In mixed hardening rule , the increments of plastic strain can be divided into two parts as 
def. = def i +def k (2.5.10) 

v v v 

Where, d£-J l is associated with the expansion of the yield surface anddef k is associated 
with the translation of yield surface. These components are assumed to be given by 
def' = Z• def , def k — (1 Z)• def (2.5.11) 

V V V V 

where Z is the parameter of mixed hardening , taking the value in the range of 0 < Z < 1 . 
It is the value Z that adjusts the contributions of the two hardening rules and that 
simulates different degrees of Bauschinger's effects. The kinematic hardening rule is 
recovered with Z = 0 while the isotropic hardening is recovered with Z = 1. 

2.7 Normalized Work in Triaxial Test 

In response to changes in the state of stress, the deformation of a body of sand can be 
divided into recoverable deformations and irrecoverable deformations (or elastic strains) 
are generally the result of deformation of the sand grains with very little or no particle 
rearrangement. Whereas, particle rearrangement is the major source of the irrecoverable 
deformations. The normalized work represents the amount of energy, per unit of mean 
pressure, required to produce a certain amount of particle rearrangement. The 
macroscopic manifestations of such a particle rearrangement are the shear strain and the 
dilatancy volumetric strain. Since particle rearrangement is basically a frictional process, 
it is therefore assumed that the normalized work should be independent of bulk 
properties, and test conditions such as the stress path and mean pressure. The normalized 
work is treated as a material property which enables the determination of dilatancy 
volumetric strain under arbitrary stress conditions and thus ideally suited for 
incorporation into general constitutive models for granular materials. 

In a triaxial compression test (cJj > 02 = 0 3 ), the stress difference q = aj — <73 and 
shear strain y = £y - £3 are convenient measures of deviatoric stress and deviatoric strain. 
Using these variables, the expression for normalized work can be written as follows. 

dSl{y p )=deP d +lldr p (2.7.1) 

3 p 
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In which, d y p = deF - deF , and p = — (<Tj + 2a ^). All the experimental results from 
1 J 3 

triaxial tests which will be presented in later in fig. 2.7.1. It can be seen that the slope of 


such a curve is 


equal to 


CO. 



Shear Strain, y 



Fig. 2.7.1 Schematic representation of normalized work for triaxial compression 

The normalized work computed from Equ. 2.7.1, will be dependent on the stress path. For 
example, the normalized work computed from the results of triaxial extension test will 
fall below that of triaxial compression test. The normalized work computed from the 
results of tests with any other stress path (for example pure shear tests) will fall between 
those of triaxial compression and triaxial extension. Due to the dominant effect of the 
second term on the right hand side of Equ. 2.7.1, it can be shown that the normalized 
work for various stress paths are approximately proportional to the shear strengths 
(ultimate value of q/p) along those stress paths. 

The variation of shear strength for different stress paths can be represented by a function 
R( 0 ) in terms of the Lode angle, 6 , defined by the following equ. 



(2.7.2) 



Where, J 


The function R( 0 ) can be incorporated in the normalized work to approximately account 
for the effect of the stress path. It is postulated that, if Q. denotes the normalized work 
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for a triaxial compression stress path, then R(#)f2 is the normalized work for any other 
stress path. As a consequence, generalized forms are; 


d£l = R(e)cQ(£)d% 


(2.7.3 ) 


dF p ~ 
at Vd 


R{0M£)--S u Eij 


P 


d% 


( 2.7.4 ) 


In these equations Q){<%) is the slope of the normalized work for triaxial compression test 
condition. 


2.8 Stress-Dilatancy Relation 

The stress-dilatancy relations, as the name implies, refer to theories which attempt to 
explain and / or model the volumetric strains during the shearing action. In a sense, the 
stress-dilatancy relations can be considered a form of plastic flow rule. Many stress- 
dilatancy equations have been proposed in the literature. Among the earliest stress 
dilatancy equations is the one proposed by Taylor (1948) for direct shear, which was 
subsequently extended by Bishop (1954) for triaxial compression. Other stress dilatancy 
equations have been proposed by, among others, Rowe (1962), Rowe, et al. (1964), 
Horne (1965), Barden (1969) and more recently by Dejong (1976), Moroto (1976) and 
Tatsuoka (1976). Most of these studies, with the possible exception of Moroto's work, are 
limited to specific conditions such as plane strain or triaxial state and the main differences 
arise from the selection of stress components or a measure of the state of stress to 
describe the stress -dilatancy behavior. 

The relationship that needs to be specified is stress dilatancy relation for the complete 
characterization of the behavior of sand. While stress strain relation tackles shear of 
distortional deformation, volumetric deformations are dealt with stress dilatancy relation. 
These kinds deformation are linked with each other and are found to depend on both 
initial void ratio and stress condition, Cubrinoviski an Ishihara, 1998. Such a link can be 
established by a relation based on energy approach, Roscoe et. al. 1963. Gutierrez, later, 
modified the original cam-clay flow rule by introducing a term ' non-coaxiality' that exists 
between principal stress and plastic strain increment direction. Dafalias( 1999 ) has given 
the stress-dilatancy relationship for the case of triaxial shearing as, 

rj = M-jU-d ( 2.8.1 ) 
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Where, d is dilatancy = 


df df P q 

v ~ v , jj = — , d£ p is plastic ( volumetric ) shear 


d ^q del! 


strain, def! is plastic ( deviatoric ) shear strain, , and jj = ^ , is rate of energy, W is 


d£ p q 


dissipated during the shear distortion of a unit volume of soil, under mean normal stress, 
P (o = ^) .Based on experiment, Ishihara and Kabilamany ( 1990 ) found that 


de p q 


increases with increase in d £ p , and modeled as; 

H 


M = Mo + — \M ~M 0 \tan 1 


d£ p 


(2.7.2) 


Where, // 0 and M are slopes of the normalized shear work verses plastic shear strain 
relationship at small and large strains respectively. S c is shear strain at which 

H - ^ is attained.To realize the steady state of deformation at shear strains, 


M = 


W' 


, cs is critical state stress ratio. 


Based on the experimental evidence of Q. - d£ p relations for sand ( Momen and 

Ghaboussis, 1982; Cubrinivski and Ishihara, 1998 ), it is assumed that the parameter of 
the stress-dilatancy relation, // 0 .M. S c are independent of both initial density and normal 
stress of sand. Effects of the initial e-p’ state on the stress-dilatancy relation are implicitly 
incorporated through the parameter, //. 

The dilatancy parameter that specifies the rate at which energy is dissipated during shear 
distortion of the unit volume of soil is given by the slope of normalized shear (Q) versus 

plastic shear strain , [e : p ) plot, fig. 2.7.1. Study on such relationship for various sands 

suggests that Q-(e^) plot is not affected by factors such as relative density, mean 

pressure, overconsolidation ratio and inherent anisotropy Momen and Ghabousi 1982. 
Experimental data leading to the above conclusion were obtained in triaxial tests on six 

different sands. This leads to the conclusion that a nearly unique Q- [e p ) relation exist for 
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all initial states of a given sand. Further, the dilatancy parameter appearing in the relation 
based on energy approach is found to increase with increasing plastic shear strain. 
Kabilamany and ishihra (1990) modeled this tendency with the help of multiple shear 
mechanism. 


2.9 Review of Undrained Behavioral Sand Model of Cubrinovski and Ishihara 

The plastics behaviour of sand is specified by a couple of basic relation: a modified 
hyperbolic stress-strain relation and an energy based stress-dilatancy relation. Both 
relations include basic parameters that are dependent on the magnitude of the plastic 
shear strain. Cubrinovski and Ishihara, ( 1998) has shown that the behaviour of sand 
under the combined influence of initial void ratio and confining pressure can be modeled 
satisfactorily using state concept over a wide range of density and pressure that are 
practically relevant. 

■ Stress-Dilatancy Relation 

The stress-dilatancy relation is based on the energy approach ( Roscoe et al 1963 ) is 
given by : 

= (2.9.1) 

deP P 

Where, q is shear stress, c is non-coaxiality term. 

■ Incremental Formulation ( Plasticity Formulation) 

Two dimensional formulation : 

Stress Increment Vectors = {(>Y = {<7 X a y T xy } 

Strain Increment Vectors = {def = {de x df v d y } 

Model formulation in X-Y-p’ plane : 
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Fig 2.9.1 Stress space and stress components 
The variables are 


X=^((Ty-(T x ) 

(2.9.2) 

II 

3 

(2.9.3) 

p ^ (cr v + <7 x ) 

(2.9.4) 

q = yl X~ +Y 2 

(2.9.5) 

q 

v=- 

(2.9.6) 

P 


d£y = d£% + d£(’ 

(2.9.7) 

d£q = -\l (d£j +d£^ ) + (d y%y )~ 

(2.9.8) 


The general constitutive equation for the elastic-plastic model, is that the total strain 
increment, de- is supposed to be sum of the elastic part, de-j and the plastic part, defj 


d£ u =d£-)+d£j] (2.9.9) 

l J v LJ 

The elastic behaviour is supposed to be linear and isotropic, and the plastic behaviour is 
anisotropic. 

The incremental elastic relations are given by 

da V] E ljkr d£' kl (2.9.10) 

Where, Ejjpj is the elasticity matrix. 
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The incremental plastic relations are given by 

d 4=( L )p- (2.9.12) 

Where, direction of plastic strain increment is specified as the normal to the plastic 
potential, g whereas the magnitude of plastic strain increment is given by the loading, L. 


■ Plastic Deformation and Failure 

It is assumed that a purely elastic region dose not exit or sand yields continuously from 
the very beginning of load application. The failure surface is assumed to have a circular 
cross section in the X-Y plane; for an isotropic material, 



Fig 2.9. 2 Failure Surface 

f = {x 2 + Y 2 Y 2 -ri max -p = 0 (2.9.14) 

f \ 

Where, // niax = q 

^ P ' max 

■ Flow Rule 

The model adopts the plastic potential formulation proposed by Gutierrez et al (1993 ). 
This formulation is derived upon particular considerations of stress paths including 
rotation of principal stresses and is characterized by the dependence of the direction of 
plastic strain increment on the direction of the stress increment. 
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Fig. 2.9.3 Flow Rule 


cl£ l y - de? - L 


dx 


X=X C Y=Y C 


L.X C 
7max P 


dT? y =L 


And, 


dY 


L.Yr 


J X=X C Y=Y C 


7max P 


(2.9.15 ) 


(2.9.16) 


f=k 2 + r c 2 p 


(2.9.17) 


Vmux ' P 0 
The straight line passing the points A and B is 

Y c = (X C -X).tan2g + Y (2.9.18) 

Where, 2q is the angle between the direction of stress increment and X-axis. 

d£% = — \p-T].c) -} 

x 7ma x-P } 


del - 
y 2 

dTl, = L 


= -{jU-TJ.c)+ Xc } 
1 ' ' Vmax-P 


(2.9.19) 


XV 7max -P 


■ Loading Surface and Plastic Modulus 

For a specific loading surface, 1, the value of the plastic modulus, H p is constant and its 

magnitude is related to the size of loading surface. For the cross section in normalized 
X Y 

— — — stress space, and an isotropic evolution of the loading surface, in which the plastic 

shear stain is used as a hardening parameter, or that control the size of surface fig.2.9.4. 
The loading surface is, 
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l = l \pij^q^=9~fl-P = 0 


(2.9.20) 


f Y/p 



Fig 2.9.4 Loading Surface 

9/ 9/ 

p= “Va< 

Simplified as, 

2.10 State Of Stress Conventional Tri-Axial Shearing Mode Tests 


(2.9.21 ) 
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2.11 Physical Parameters of the Bagmati Sand 


e max 

e min 


= 0.946 
= 0.610 


D 50 = 0.6 mm 
Cu =2.65 


C C -112 

% Fines = less than 1 % 
G s =2.68 


Source: Jha,2004 


Chapter 3 METHODOLOGY 
3.1 Stress-Strain Representation 

A major simplification applied to the constitutive to be developed is the limitation of the 
formulation to two dimensions. The cylindrical specimen of the conventional triaxial test 
is in axisymmetric stress state. It is considered that the assumption of a two dimensional 
condition will avoid complicating the formulation significantly and yet on the other hand, 
will be appropriate for idealization of many in-situ loading condition within reasonable 
computational efforts and required accuracy. 
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Fig.3.1.2 Representation of monotonic triaxial compressional on II -Plane 

For the shearing mode, o 2 = c 3 and r 13 , and the stress point moves along the 

(<j ] - (J 3 )- axis. When it touches f m , for the conventional triaxial compression condition, 

cr, - 0-3 = a™ + C (3.1.1) 

Where, a" 1 and k"‘ are the initial position and size of the yield surfaces, respectively. 

3.2 Incremental Formulation 

The model is formulated in the X-Y-p space, where the stress variables are defined as; 


X = 


V3 


(CT| +2-CT 3 ) 


(3.2.1) 


42 














(3.2.2) 

( 3.2.3 ) 

(3.2.4) 

(3.2.5) 



P 

(a) (b) 


Fig.3.2.1 ( a ) Stress Components , ( b ) Failure surface ( K,H, Roscoe, A.N. Schofield, 
1961 ) 

The incremental vertical strain is given by, 



(3.2.8) 


Where, G e and FIp are Elastic Shear Modulus and Plastic Modulus respectively. 
For the convention triaxial test with axisymmetric stress state, 
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(3.2.9) 


2 

de q =[de n \--\de n -de 33 ] 

Since , in case of undrained shearing , de v is zero, which suggests de q = d£\ ] 

i.e. axial strain is similar to the shear strain. 

To start, it is assumed that the total strain increment d£{j is decomposed into elastic and 

plastic components d£f- and d£- 3 - by a simple superposition , 

L J l J 


d£jj = d£jj + dejj (3.2.10) 

The elastic and plastic strain increment components are calculated separately, by an 
elastic stress-strain law of the form, 


rfcTy — Djjfci ■ d£^j ( 3.2.11 ) 

Where, 

Eijkl = | _ &ij ' dkl + (’e i^ik ' $jl + il ' & jk ) ( 3.2.12 ) 

v is Poisson’s Ratio, S(j = Kronecker Delta (Sjj =0 if i & j , and = 1 for i = j ),G e is 

Elastic Shear Modulus, dependence on the e-p-plane and empirically defined as (Iwasaki 
and Tatsuoka, 1977 ) 


Ge — 2E Pa 


(2.17-g) 2 
1 + e 


P 

Pa 


n 


(3.2.13) 


Where, A is dimensionless parameter, p a is atmospheric pressure (lx lO 4 Kpaj, p is the 
mean normal stress, n is the dimensionless material parameter determining the variation 
of G e with p. 


The incremental plastic relations are given by 

d£fj (3.2.14) 
3 d<J;: 


Where, direction of plastic strain increment is specified as the normal to the plastic 
potential, g whereas the magnitude of plastic strain increment is given by the loading, L. 
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(3.2.15) 


L 



Kdp+V.Jg 

dp dq 


p dp 

H P 


And, 


H P = 


G N~ f ■ przr ( G N - G N, min ) 


0.01 



[ r ] 

>x- 

1 _ 

1 


r fl 


X P 


(3.2.16) 


Where, Hp is plastic modulus, and 1 is loading surface. Along the loading surface, plastic 
modulus is constant. The loading surface which is moving in the stress space, a(j, 

according to the load increments, is used to evaluate the current magnitude of the plastic 
modulus as well as to trace the loading intensity. It has to notice that the G N is the initial 
modulus whereas, Hp is the tangential plastic modulus. 

Finally, the general expression for the plastic strain increment can be written as, 


deE= ‘ 


l J H P 


3a,j V 


dg 

da;; 


The total strain increment can be calculated as 


d£ij - D 1 ijkl ■ da^j + 


H P 


£ij da ' li 


_s$_ 

da;; 


(3.2.17) 


(3.2.18) 


Since, our experiment is strain controlled, so its necessary to model a strain controlled 
simulation. 


da ij= D ijkl\‘>e kl 


Hp 


— da. 


dg 

da;; 


(3.2.19) 


The calculation of the plastic strain increment components in the framework of the 
incremental theory of plasticity is based on definition of the yield function, plastic 
potential function, and hardening function. The stress-dilatancy relationship for the case 
of triaxial shearing, given by Dafalias( 1999 ) is used in the formulation. 
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Fig. 3.2.2 Representation of Stress Space in Triaxial Compression Loading. 

In the principal stress space, the yield surface is, 

0 (3.2.20) 

And, 

The failure criterion is 

/ = | • (sy - of )■ (s /7 - of)+ r 2 ax = 0 ( 3.2.21 ) 

Where, S it ( = - p ■ S n ) is devistoric stress tensor, CC'" is the coordinate of the center 

v L J U l J 

of the yield surface, f m . 


% -C/ -7max 

V = --o ’ 

%iax ^/V -£q 

^max — a l ^1 Ay > As — ^ 

Where, ai and bi are stress strain material parameters. 


(3.1.22) 
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3.3 Material Parameters Determination 


A distinct feature of the constitutive model of the present study is that a single set of 
values of the material parameter is used to model the behaviour of the given sand. That is 
whatever is the relative densities of the sand and the confining stress, the values of the 
material parameters remain the same. Thus, it is a prerequisite to use a determination 
procedure and to evaluate material parameters which will able resonablely accurately 
overall representation of a given sand. In general, except for the elastic parameters, the 
model parameters are dependent on the fabric of packing. And, except the mean normal 
pressure of the reference lines, all the parameters are dimensionless. The model requires 
four groups of parameters. 

a. Elastic parameters. 

b. State lines (reference states). 

c. Stress strain parameters. 

d. Dilatancy parameters. 

a. Elastic parameters. 

There is no specific requirement for the evaluation of the elastic parameters, and 
insistence on determination of these parameters at very small strain levels is not necessary 
since modeling of small strain behaviour is not attempted herein. Moreover, the Poisson’s 
ratio and the exponent, n may be assigned fixed value of 0.2 -0.3 and 0.5 -0.6 
respectively. Thus reducing the problem of evaluation of parameters which can be back 
calculated from known value of elastic shear modulus, G e ; equation 2.4.13. 


The empirical relation for G e in undrained case is given by (Hicher, Pierre-Yves, 1994 ); 


G e = 


_ 450V7 (1 + vJ 


1 + v 


(3.3.1) 


Where, v u is the Poisson’s ratio in undrained condition (=0.5) for perfectly saturated soil. 


b. Reference states. 

The state lines, the UR-lines and QSS-line, will be evaluated using results of several 
monotonic undrained strain controlled tests on loose samples having different void ratio 
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with the only requirement being that in the course of the undrained shearing, the samples 
will exhibit drop in the shear stress or QSS. Assembling the quasi steady states from each 

test in the e-p’ diagram, QSS- line will be drawn. If during the test of undisturbed soil 
samples, the requirement for drop in the shear stress often may not be satisfied, and in 
such case it will be necessary to replace the QSS-line with SS-line as a reference line. 

Again, if the zero steady state can not be directly observed in the undrained tests, the 
position of the UR-line (threshold void ratio ) will be reasonably approximated in the e- 
p’ plane by extending either the QSS/SS -line to the point of the intersection with the 
void ratio axis. This point of intersection then denotes the void ratio ejj. 

c. Stress-strain parameters 

The stress strain parameters of the model consists of linear coefficients of the expressions, 

Vimix — Cl\ + b\I s , Is — 0 

E N, max = a 2+ 0 ( 3.3.2 ) 

E N, min ~ a 3 + ^3As > — 0 

The coefficients, a\,b\,a 2 ,b 2 ,a^,b^ will be determined using results from a series of 
drained p-constant tests on samples with distinct relative initial states or samples having 
different values of the state index. 

The peak stress ratio, // max of each test will directly read off from the measured 
normalized stress strain curve and the initial plastic moduli, E N max and E v min will be 

determined either as the best fit small strain and large strain values respectively, in trial 
and error simulations of the experimental data or by establishing an experimental 

relationship between E N and d£[] ( Cubrinovski and Ishihara, 1998 ). For each test, the 
state index I s will be calculated using the reference line and the initial e-p state of the 
sample, and the values of ?7 max , E Wmax and E Vmin as defined above will be plotted verses 
I s ; these correlation with I s will be approximated with linear relationship and the 
coefficients, ci\ , b\ . a 2 ,1'h , c /3 , £>3 will be determined from these data arrangements. 
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d. Dilatancy parameters 


For the evaluation of dilatancy parameters, the critical state stress ratio, M =« will be 

l p) cs 

first determined from the aforementioned monotonic undrained tests or alternately it may 
be approximated by the peak stress ratio of the drained p-constant tests, for I s = 1, or the 
final slope or slope at large strains of the curve, normalized plastic shear work versus 
plastic shear strain from drained tests. The parameter jj. 0 is the initial slope or slope at 
very small strains of the curve, normalized plastic shear work versus plastic shear strain 
from drained tests. The dilatancy parameter, S c which controls the pore pressure 
development and has an especially large influence on the cyclic strength. The parameter 
will be determined by mathematical corroboration due to absence of cyclic triaxial 
equipment in the lab, or can be determined as the value of the plastic strain at which 
average value of // is attained. For a larger value of S c reflects smaller /i at a given 
strain or less contractive behavior. The material parameters for the Bagmati sand will be 
used to model the monotonic behaviour of samples with various combinations of initial 
relative densities and confining pressure. 
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Chapter 4 RESULTS AND DISCUSSION 
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APPENDIX - A 


Computation of Normalized Work 

In general, the normalized work can be computed from the results of any experiment 
which induces deviatoric stresses in the sample, including standard triaxial test. Both 
stress and strain measurements are required. Moreover, if the mean pressure during the 
experiment is not constant, the results of an isotropic compression test are also required in 
order that consolidation volumetric strains can be determined. 

If the elastic material properties are available, or can be easily determined, then the elastic 
components of strains can be subtracted from the total strains to determine the plastic 
components of the strains. Using these plastic strains, the normalized plastic work can be 
computed. However, the elastic components of strains are usually a small part of the total 
strains. It is often convenient to use the total strains rather than the plastic strains. In this 
way the normalized work is computed and usually it is a very close approximation of the 
normalized plastic work. The differences between normalized work and normalized 
plastic work are often negligibly small, for all practical purposes. All the experimental 
results to be presented in later sections are in the form of normalized work, computed 
from total strains. 

In order to calculate the normalized work from the results of a standard triaxial test, first a 
plot of the stress ratio q/p versus shear strain, y = £ \ — £ ^ > is prepared. Then the 

following steps are performed. 


1. Divide the shear strain axis into a number of increments, Ay,and then for each 


shear strain increment calculate the area under the curve [ ~ , A y ]. 


2. Determine the increment of volumetric strain, Ae v — A.£\ + 2A £'3 

3. Determine the increment of consolidation volumetric strain. 

. _ C c A p 

A£VC - 2.3(1 + e 0 )P 

<?0 = initial void ratio, C c - Coefficient of compressibility, determined from isotropic 
compression test. 
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In a constant cell pressure test A /7 =A(Xi/ 3 ; otherwise, in general 

A _(Acr 1 + 2Acr 3 )/ 

P ~ /3 ‘ 

4. Determine the increment of the dilatancy volumetric strain. 

v — vc 

5. The increment of the normalized work is computed from the following equation. 

A£l = Ae vd +=--Ay . 

3 p 

The increments of normalized work are summed and plotted versus the shear strain, y. 
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APPENDIX - B 


Computation of Elastic and Plastic Strains 
Elastic Part 

For a given elastic strain increment, the stress increment can be calculated as, 

{cl of = [£]{/(/'f 


For plane strain condition elasticity matrix is given by 


[E] 


Ej 1-v) 

(l + v) • (l - 2v) 


1 


v 

(l^v) 

1 


(1-v) 

0 0 


0 


0 

(l~2v) 

2-(1-v) 


A1 

A2 

0 

A2 

A1 

0 

0 

0 

A3 


Where, 


.u- e-d-b 

(l + v).(l-2v) 


A2 = 


E-v 

(l + v)(l - 2v) 


A3 = 2 G e 


Hence, expression for stress increment in simplified form can be expressed as 
do ^ = A1 • de I + A2 ■ de* 

doy = A2 ■ de 3 + A1 • de{ 

d o 12 = A3 ■ de{^ 


The elastic shear constant , G e , can be calculated from the equation (Iwaski and 
Tstsuoka, 1977) 


G e = A- p a 


(2.17 - e)“ 

(l + e) 


\ n 


Pa J 


Where A, n are known as elastic shear constant and exponent that determines the rate of 
variation of G e with p respectively . In this equation , pressure of mean normal stress and 
void ratio terms signify that the effect of initial state is also included in the elastic part. 
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Plastic Part. 


To derive the expression for incremental plastic strain, the expression for loading index 
can be expressed in matrix form as 

Y 

[L] = — [Bl B2 B3][da] T 

H P 


Equivalently, 

tip 

Where, 


B\ — — 
2 




B 2 = - r + 


X r 


r f 


B3 


7 


v j y 

Again, 

[daf = [E]{de} T 

= \E]{le-d£ p J 

= \E\de} T -(d)zJ^ 



d(7 u 


And, 


L = ^-{B\[Etde} T -[ e\l\ 


_9g_ 

doij 


On Simplification gives 
L= {BjFjde} 

H ' +{BW fe} 
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Where, 


[ ^ ] 



7 M -7^ 

Xc 


d (Tj 



{ B ) 

7 


dg 

, - j 

1 

-Jj] 

Ac 
_|_ 


da 3 

> — 5 

2 " 

{ B ) 

7 


dg 



Yc 



_ 3<7 13 J 



7 

J 



Cl 

C2 

C3 


1 


A little consideration will show that 


{B\E\cle} T =Klde 3 +K2d£ l +K3de l3 
Similarly 

-p: j = KICl + K2C2 + K3C3 

Where, 


Al = A151 + A252 
K2 = A\B2 + A2Bl 
K3 = A3B3 


Finally, the expression L can be written as 

K 1 • d£^ + K2 ■ d£ 3 • +K3 • d£^ 3 
H p +K1C1 + K2C2 + K3C3 

This expression for L is suitable particularly when stress increment is calculated from 
strain increment. Now for a given total strain increment, the plastic strain increment can 
be calculated as. 

d£% -L C\ 
d£[ -L C2 
d£[ 3 =LC3 
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APPENDIX - C 



STOP : If Negligible Stress Increment 
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